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(54) Dynamic light scattering method and apparatus 



(57) A dynamic light scattering apparatus comprises 
at least two lasers (11) of different wavelengths which 
illuminate the same volume (SV) of a fluid containing 
scattering or fluorescent centres. The scattered or fluo- 
resced radiation is detected at each discrete wavelength 
and properties, such as the concentration of different 
particle sizes, of the volume (SV) are derived. The de- 
tection at multiple wavelengths allows the effects of mul- 
tiple scattering to be overcome. Multiple volumes (SV) 
may be illuminated and detected simultaneously so as 



to measure the spatial variation of particle size concen- 
tration. The rotational diffusion constants of the particles 
may be measured from the depolarized component of 
scattered light. 

The apparatus or method may be used to measure 
the velocity of particles in the volume (SV). The appa- 
ratus or method may also be applied to remote sensing 
of the atmosphere by scattering of laser pulses, the in- 
tensity variation functions of the received scattered puls- 
es being concatenated prior to auto- or cross-correlation 
to detect particle size concentrations. 
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Description 

[0001] The present invention relates to a dynamic light scattering (DLS) or photon correlation spectroscopy (PCS) 
method and apparatus, particularly but not exclusively for Fluorescence Correlation Spectroscopy (FCS) or lidar ap- 
plications. 

[0002] PCS/DLS involve the illumination of a sample with coherent light. The intensity of photons scattered from the 
sample fluctuates as a result of the Brownian motion of scattering particles in the sample. The number of photons 
scattered in a narrow range of angles, defined by specific apertures and light propagating and focusing elements, is 
detected as a function of time and the resultant function is autocorrelated, preferably in real time. This allows the 
distribution of relaxation time constants and therefore the distribution of sizes of the scattering particles to be deter- 
mined. 

[0003] FCS involves the use of a similar technique on fluorescent samples, except that the photons emitted by flu- 
orescence, instead of scattered photons, are detected as a function of time. The resultant function is autocorrelated 
in real time, in order to measure the specific change in binding of molecules, detect molecules present in low concen- 
trations and determine related reaction rates. 

[0004] DLS-Lidar involves the active remote sensing of a portion of the atmosphere by illuminating the portion with 
coherent light and detecting the intensity of photons scattered from that portion as a function of time. The resultant 
function is autocorrelated in real time, in order to derive the particle size distribution within the sampled portion of the 
atmosphere. 

[0005] A problem associated with the above techniques is that, as the number per unit volume of scattering particles 
increases, there is an increased probability of multiple scattering in which a detected photon has been scattered by 
more than one particle, which means that the desired properties can no longer be measured directly. Furthermore, the 
range in velocity of Brownian motion becomes restricted because of collisions between scattering particles. In order 
to compensate for these effects, the normal practice is to detect the scattered photon intensity at a range of different 
scattering angles, typically from 10° to 180° (backscattering). As there is a quadratic dependence of time constants 
on scattering angle for the ideal case, any deviation from this quadratic function, due to multiple scattering or restricted 
particle collision, is quantifiable and can be compensated for. However, such procedures are time-consuming because 
of the number of readings that need to be taken, and require expensive goniometers. 

[0006] In lidar applications, it is not practically feasible to perform angle-dependent scans of a portion of the atmos- 
phere, so that the particle size distribution and degree of multiple scattering cannot readily be measured. 
[0007] According to one aspect of the present invention, there is provided a dynamic light scattering apparatus and 
method in which a sample is illuminated by laser light of different wavelengths and the fluctuations in light scattered 
or fluoresced by the sample at each wavelength are detected. The time constants of particles in the sample are derived 
from the detected light fluctuations by auto-correlation or cross-correlation of the detected light intensity with respect 
to time. Measurement of light scattering or fluorescence at different wavelengths replaces measurement at different 
scattering angles, thus simplifying the process of measurement. 

[0008] According to another aspect of the present invention, there is provided a method and apparatus for measuring 
the velocity of particles in a sample by means of dynamic light scattering, in which the frequency of oscillations in 
correlation functions of the detected intensity fluctuations as a function of time are measured and used to determine 
components of velocity in the direction of detection. 

[0009] According to another aspect of the present invention there is provided a method and apparatus for the meas- 
urement of properties of particles in a section of atmosphere by means of dynamic light scattering, in which the section 
is illuminated with laser pulses and the pulses scattered from that section are detected and their intensity fluctuation 
correlated to determine the properties of the particles. Preferably, the intensity fluctuations of the detected pulses and 
concatenated and/or high-pass filtered prior to correlation. In this way, dynamic light scattering techniques may be 
applied to remote sensing of the atmosphere. 

[001 0] In an embodiment of the invention, a dynamic light scattering apparatus comprises discrete numbers of trans- 
mitting lasers with various wavelengths ranging from the visible to the near infrared. Each discrete laser transmitter 
forms a transmitter channel, which operates in conjunction with one or more receiver channels. The laser light propa- 
gates to the light scattering sample via apertures, focusing, telescope and/or fibre optical elements. By these elements, 
beam waists are produced in the sample. Corresponding to the channels of transmitting lasers are a second set of 
functionally identical elements, for collecting and receiving each wavelength of scattered light separately and simulta- 
neously. The apertures, light collecting and propagating elements for one channel are adjusted in combination to pro- 
duce an intercept of the dynamic part of the autocorrelation function close to unity. 

[0011] The automated functions of the apparatus from one back scattering position replace the function of a goni- 
ometer for scattering angle scans. This is achieved by (quasi-) simultaneous, n-fold channel operation. The apparatus 
enables classification and quantification of multiple scattering sample systems, the normalisation of particle cross- 
sections, sizes and shapes to the wavelengths used and the quantification of repulsive forces exerted to particles or 
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restricted particle collision as in a gel-like system with cage functions. If multiple scattering is encountered the apparatus 
allows adjustments to be selected to derive the true particle size. The simultaneous multiple volume, n-fold, acquisition 
capability of the apparatus applies to three general fields: particle sizing in the atmosphere by active remote sensing, 
particle sizing and visco-elastic property determination in fluids and suspensions as for example in crystal growth of 
5 bio molecules from solutions and property determination of light emitting samples by fluorescence recording as for 
example in immune assay developments. The apparatus utilises cross-correlation, depolarised and polychromatic 
dynamic light scattering in pulsed and continuous wave operation and filtering methods to derive appropriate autocor- 
relation functions. 

[0012] Specific embodiments of the present invention will now be described with reference to the accompanying 
10 drawings, in which: 

Figure 1 is a schematic diagram of the optical front end of one channel of an apparatus used to measure back- 
scattering from suspended particles; 

Figures 1a and 1b show alternative examples of the optical front end in Figure 1; 
15 Figure 2a shows an example of recorded intensity fluctuation as measured by the apparatus of Figure 1 ; 

Figure 2b shows the autocorrelation function of the fluctuation of Figure 2a; 
Figure 2c shows the particle radius distribution derived from Fjgure 2b; 

Figure 3 shows the variation of time constants with scattering angle for different particle sizes in a sample of 
aqueous solution or fluid suspension with a given viscosity and refractive index; 
20 Figure 4 shows the variation of time constants with scattering angle for different particle sizes in a sample of air, 

gas or atmosphere with a given viscosity and refractive index; 

Figure 5 shows different autocorrelation functions and bulk velocity measurements from particles in samples of 
fluid at different velocities; 

Figure 6 shows autocorrelation functions and corresponding radius distributions from samples with contained par- 
25 tides having a relative velocity to each other, showing oscillations of autocorrelation functions prior to nucleation 

of these particles to form a crystal; 

Figure 7 shows the measured angle-dependent deviation of decay time constants from a theoretical curve as 
shown in Figures 3 and 4 with the given viscosity and refractive index; 

Figure 8 shows the variation of scattering volume length with atmospheric height at different separations between 
30 transmitter and receiver; 

Figure 9 shows the equivalence of time constants and wavelengths with scattering angle as calculated and nor- 
malised for two different wavelengths; 

Figure 10 shows a generalised embodiment of the present invention for the measurement of a sample; 
Figure 10a shows beam tracings in an implementation of the embodiment of Figure 10; 
35 Figure 10b shows an implementation of the embodiment of Figure 10 for measuring the depolarised cross-corre- 

lation DLS; 

Figure 11 shows the measured values of time constants for two different wavelengths; 

Figures 12a and 12b shows two alternative embodiments for active remote sensing of the atmosphere; and 

Figures 1 3a to 1 3f show the signal processing steps for measuring particle radius distribution in the embodiments 

40 of Figures 12a and 12b. 

Figure 1 shows in schematic form the optical front end of one channel (i.e. one transmitter and one receiver) of 
the apparatus consisting of a transmitting optical fibre or fibres 4 connected to a laser source and a receiving 
optical fibre or fibres 2 connected to a detector, such as a photomultiplier tube or avalanche photo-diode, mounted 
together so that the axes of their ends are parallel (as in Figure 1a) or at a defined angle (as in Figure 1b) to each 

45 other. Laser light emitted by the transmitting optical fibre(s) 4 diverges over a beam B while the receiving optical 

fibre(s) 2 collects light emitted from a reception volume RV. In a scattering volume SV which is the intersection of 
the illuminating beam B and the reception volume RV, particles scatter the laser light in different directions and 
the receiving optical fibre(s) 2 collects light scattered within a narrow range of directions. 

so [001 3] As shown in Figures 1 a and 1 b, the defined angles between the transmitting f ibre(s) 4 and the receiving fibre 
(s) 2 are adjusted by movement of tight-fitting coaxial micro-positioning stages (MPS) in the reciprocating direction 
shown by double arrows, so as to locate the scattering volume SV at distances L1 and L2 in Figures 1a and 1b re- 
spectively, and to compensate for the effect of Snell's law on different wavelengths when co-locating the beam waists 
of different wavelengths. Additional beam forming is produced by suitable optical elements (OE) such as lenses, prisms, 

55 wavelength filters and/or telescopes, according to the specific application. 

[0014] The receiving and transmitting optical fibres 2 and 4 are selected as either mono-mode or few-mode (i.e. 
more than one but preferably less than about two modes) fibres and may either receive and transmit through the same 
optical elements or have different optical elements. Different possible combinations will now be discussed with refer- 
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ence to Fig 1c. 

[0015] In arrangement I, the transmitting fibre 4 comprises a single mono-mode fibre and the receiving fibre 2 com- 
prises a single few-mode fibre (shown as option la), or the transmitting fibre 4 comprises a single few-mode fibre and 
the receiving fibre 2 comprises a single mono-mode fibre (shown as option lb); in either case, the transmitting and 

5 receiving fibres share the same optical elements. 

[0016] In arrangement II, the transmitting fibre 4 and the receiving fibre 2 have separate sets of optical elements. In 
option a (indicated by the reference ila in Figure 1c), the receiving fibre 2 is a few-mode fibre and the transmitting fibre 
4 is a mono-mode fibre, as in arrangement la, while in option b (indicated by the reference lib in Figure 1c), the trans- 
mitting fibre 4 is a few-mode fibre while the receiving fibre 2 is a mono-mode fibre, as in arrangement lb. 

w [001 7] In arrangement III, the transmitting fibre(s) 4 and the receiving fibre(s) 2 have separate sets of optical elements. 
The transmitting fibre 4 is a single mono-mode fibre, while the receiving fibres 2 comprise at least two mono-mode 
fibres in parallel, each having a separate detector, discriminator and pulse former. Preferably, the receiving fibres 2 
are arranged in pairs of fibres arranged to receive the same wavelength, with the outputs of each pair then being cross- 
correlated. This cross-correlating arrangement eliminates detector dead-times, after-pulsing and dark counts and will 

*5 be described in more detail with reference to Figure 10b. 

[001 8] The detector outputs a current proportional to the rate of photons received by a receiving fibre 2, which current 
varies with time as shown for example in Figure 2a. The current is converted to a digital value by an A/D converter 
connected to a computer, which stores a set of digital values as a function of time. Alternatively, individual received 
photons are counted by means for example of a photomultiplier tube, a discriminator for setting a threshold for the 

20 signal of a detected photon, a pulse former for generating a pulse when the photon is detected, a counter for counting 
the generated pulses within a repeating period and a storage device for storing the count values for each period. The 
computer performs an autocorrelation function on the stored values as follows: 



25 



40 



45 



50 



ACF(r) = X '(')/(' + r) (1) 



where l(t) is the fluctuating received light intensity as a function of time t, such that J(t) equals the average scattered 
30 intensity over the integration time, and x is the sampling interval, which forms quasi-logarithmic channels in the struc- 
tures of the correlator. A sample autocorrelation function is shown in Figure 2b. 
[0019] The ACF function is related to the particle sizes by the following equations: 

35 ACF(z) = g {2) = UV-e { - 2rAx) (2) 

where g( 2 ) is the second order laser light intensity fluctuation, "1 U is an arbitrary baseline, (3 is a quality factor having a 
theoretical maximum of 1, Ax is the entire sampling time structure of the correlator, and 



x c = 1/r=1/(q 2 x D) (3) 

where t c is the time constant, D is the diffusion constant and q is the scattering vector, where 

q=(47i/a)sin(0/2) (4) 

with 0 being the scattering angle and X the wavelength of the impinging light. 
[0020] The diffusion constant D relates, via the Stokes-Einstein equation, a geometrical factor, the viscosity r| and 
the hydrodynamic radius r h to the thermal energy, which is Boltzmann's constant k multiplied by the absolute temper- 
ature T as follows: 



[0021] Under ideal conditions, the time constants x c are directly related to the radius r h and to the viscosity of the 
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suspension of scattering particles. The term p (also called the intercept) is inversely proportional to the number of 
observed coherence areas (number of speckles) N, which fluctuate due to Brownian motion. The measured ACF is 
formed from many different individual functions of the form of equation (2), for each of the different particle sizes. The 
conversion of the composite ACF shown for example in Fig. 2b to the relative amplitude intensities for different particle 

5 sizes shown in Figure 2c may be performed by a suitable computer program such as is described in 'CONTIN: a general 
purpose constrained regularization program for inverting noisy linear algebraic and integral equations', Provencher, S. 
W., Computerphysics Communication 27, pp. 229-242, North Holland Publishing Company, 1982. 
[0022] As a result of the above relationships, DLS can be thought of as a point-like measurement technique with 
virtually no spatial dimension. Three dimensional (voxel) resolution requires an embodiment using a plurality of chan- 

10 nels. The spatial mode of the received light defines N together with its inverse proportionality to the number of scattering 
centres in the scattering volume. 

[0023] The configuration of the apparatus used in this embodiment includes a minimum of two mono-mode receiver 
units, in order to obtain the maximum intercept value and to enable signal processing with cross-correlation for the 
various polychromatic embodiments. The cross-correlation function is given by: 

15 

CCF(r) = £/(/)./(/ + r) (6) 

20 

where l(t) is the intensity measured by one receiver channel and J(t) is that measured by another receiver channel. 
Equation (6) is based on equation (1), but with l(t + x) replaced by J(t + x). 

[0024] Thus, the distribution of particle radius sizes in a sample may be determined from the autocorrelation function 
(ACF) , as shown in Figure 2c, or from the cross-correlation function (CCF) as in equation (6). 
25 [0025] The relation between the decay time constant, which depends on the particle size, and scattering angle is 
shown in Figure 3 for a fluid with a viscosity of 1.98 X10* 3 Nsrrr 2 (1 .98 cP) and refractive index n=1.33, and in Figure 
4 for particles in air with a viscosity of 1.71 X10~ 5 Nsrrr 2 (0.0171 cP) and a refractive index n=1.0. 
[0026] The variation of time constants as a function of scattering angle shows that backscattering measurements at 
high scattering angles are preferable to scattering measurements at low scattering angles for the following reasons: 

30 

1) A 1° error in scattering angle gives a measurement error of approximately 17.2 % at a scattering angle of 10°, 
but a measurement error of only approximately 0.2 % at 170°. 

2) Assuming that a measurement time of 1000 time constants t c is required to give a reasonably small standard 
deviation of about 0.3 %, for a 1 ujti particle a measurement time of about 5.8 minutes is required at a scattering 

35 angle of 170° , but 11.6 hours at 10°, about 120 times as long. 

3) The light intensity scattered in a forward direction increases with particles above about 200 nm in size to about 
10 times the backscattered intensity. The probability of getting unacceptably high count rates is therefore greater 
for forward scattering. 

40 [0027] Dynamic Light Scattering may also be used to measure the velocities of particles in a fluid. In one example, 
the laser light is incident on a capillary tube through which fluid flows. Light is reflected from the capillary wall onto the 
receiver and is also scattered by the fluid. The reflected light then combines with the scattered light at the receiver to 
produce interference. 

[0028] Flow velocities are derived from the time constant x c according to the following equation: 

45 

Af=1/(2*T C ) (7) 



[0029] The velocity component v in the direction of observation is given by: 

50 

= c 
v " (1-cos(9))n f 



where c = velocity of light 

55 

n = refractive index 

9 = observation angle relative to the velocity vector, and 
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f = frequency of the laser light. 

[0030] Figure 5 shows the ACPs and velocity profiles measured from a sample with flow rates of 0 to 25 mm/s. 
[0031] For more detailed velocity measurements, the velocity profile and the extent of the measured sample need 

5 to be taken into account as indicated above in the discussion of the point-like measurement character of DLS and the 
advantage of scaling up systems to a multi-channel (n-fold) approach enabling the resolution of phenomena in three 
dimensions. The interference between the reflected and scattered light (heterodyning) can be used to resolve velocity 
gradients, but extracting the velocity of each of the n-fold channels simultaneously provides a complete profile acqui- 
sition on the basis of homodyning. 

w [0032] Pronounced oscillations of ACF's are detected prior to nucleation as a result of currents or velocity gradients 
within a suspension. Figure 6 shows the oscillating ACF (ACF-33) close to the baseline. The time constant of oscillation 
is 185 ms giving a frequency of approximately 5 Hz which is used in the equation (8) as Af, so that the velocity gradient 
in the volume can be determined. This oscillation is no longer present in a measurement of an ACF taken approximately 
10 minutes later (ACF-40). 

75 [0033] Under ideal DLS, with no particle interaction, light is collected from one coherence area (1 x N) of the scat- 
tering volume and is analysed by autocorrelation. Using two detectors receiving scattered light from the same area 
reduces the effect of noise and imperfections of the detector by cross-correlating the signal from both detectors, and 
multiple scattering is suppressed if one coherence area is observed. 

[0034] The effect of multiple scattering and deviation from free diffusivity is determined conventionally by an angle- 
20 dependent test on a sample having a known distribution of particle sizes, using a goniometer For example, Figure 7 
shows the measured deviation from the theoretical curve (dotted) of a sample with a particle radius of 392 nm, with 
viscosity of 9.0 x 10" 4 Nsnr 2 (0.90 cP). 

[0035] If a measured sample includes fluorescent particles, the sample may be illuminated simultaneously with light 
of a number of discrete wavelengths and light corresponding to the wavelength of fluorescence is detected as described 
25 in the third mode-selecting arrangement described above. The detector includes a polychroic element to separate out 
the desired wavelength. 

[0036] In lidar remote-sensing applications, a portion of the atmosphere is illuminated and backscattering from that 
portion is detected. The scattering volume is defined by the intersection between the field of view of the transmitter 
and the receiver. The scattering volume length (svl) varies as a function of atmospheric height and separation a between 
30 the transmitter and the receiver. An example of this variation with a 0.1 mrad full angle receiver aperture is shown in 
Figure 8. 

[0037] As it is very difficult to observe the same portion of the atmosphere at different viewing angles, lidars are 
usually restricted to a single backscattering angle, which restricts the number of properties which can be measured. 
With pulsed illumination, backscattered light can be resolved according to time of flight. The temperature profiles in 

35 the atmosphere can be measured from elastic backscattering signals and from Raman effect return signals from ni- 
trogen. The air viscosity can then be determined from the air temperature at the relevant height. 
[0038] In embodiments of the present invention, in order to overcome the problems associated with angle-dependent 
measurements, measurements are made instead at multiple different wavelengths while the transmitter and the re- 
ceiver remain in the same position. The deviation of the measured from the theoretical time constants i c is measured 

40 at each of the wavelengths, so as to determine the effect of multiple scattering or the degree of particle interaction 
which leads to a decrease in free diffusivity. Figure 9 shows the equivalence of time constants t c and wavelengths as 
calculated and normalised for two different wavelengths of 266 nm and 690 nm. The numbers shown next to the points 
of the graph correspond to the equivalent wavelengths, in pm, that would produce the same time constant t c as the 
scattering angle shown on the abscissa and the actual wavelength (266 nm or 690 nm), if the equivalent wavelength 

45 were to illuminate the sample at an angle of 180° (backscattering). The scattering particle size is 2 nm. 

[0039] Different embodiments of the present invention will now be described for different applications of the present 
invention. 

Particle Size and Multiple Scattering Measurement 

50 

[0040] Figure 10 shows an apparatus according to an embodiment of the invention, for the measurement of particle 
size, particle collision and multiple scattering effects in a sample. There are a plurality not transmitters 10 and corre- 
sponding receivers 12. By means of appropriate optical elements (OE), as shown in Figures 1a and 1b, the scattering 
volumes of each transmitter/receiver channel may be superimposed so that the measurements of each channel are 
55 made on the same scattering volume, or the scattering volumes of the channels may be arranged in an array, so as 
to provide two or three dimensional resolution of the measured properties. 

[0041] Each transmitter 10 comprises one or more lasers 11 each coupled to a respective mono-mode fibre 4, for 
transmitting a single wavelength X v A^, A3... Each receiver 1 2 comprises a plurality of mono-mode fibres, each having 
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a corresponding detector and discriminator, as in arrangement III described with reference to Figure 1c. The mono- 
mode fibres may be selected to transmit respective different wavelengths or some of the fibres may transmit the same 
wavelength. 

[0042] The detectors for each receiver 1 2 are shown collectively as reference 1 4 in Figure 1 0, while the discriminators 

5 are shown collectively by reference 16 and the pulse formers by reference 18. For each channel 1 to n, there is a 
correlator 20a to 20n which auto-correlates the signal from one detector, or cross-correlates the signals from two or 
more different detectors of light of the same wavelength. The outputs of the cross-correlators, in the form of cross- 
correlation functions CCF as derived in equation (6) above are input to a processor 22. In a typical implementation, 
the discriminators 16, pulse formers 18 and correlators 20 are provided by a dedicated card installed in a general 

10 purpose computer which is represented by the processor 30. 

[0043] The apparatus shown in Figure 10 allows measurements of the time constant t c at multiple wavelengths 
simultaneously from the same sample. In the example shown in Figure 11, measurements of the time constant t c for 
particles of radius 15 nm are shown for wavelengths of 632.8 nm and 689.5 nm, with their equivalent scattering angles 
normalised and calculated for a beam of 61 8 nm impinging at 1 80°. Four measurements are made for each wavelength 

15 and the maximum, minimum and average measurements being shown. The large systematic error for measurements 
at 632.8 nm are due to temperature control instabilities in the laser diode used for that wavelength, leading to limited 
coherence lengths or spectral purities. The measurements are shown in Figure 11 on a graph of time constant t c 
against scattering angle to show that the equivalencies exist and the sample has the ideal properties within the standard 
deviation for these tests. The equivalence is calculated by combining equations (3) and (4) and normalising this to an 

20 assumed wavelength impinging at 180°. 

[0044] The array arrangement allows for particle sizing as a function of extended length in a cuvette or a micro-vial 
plate. To ensure test reproducibility and the homogeneity of the sample, control of sample properties is provided by 
the quasi-simultaneous sampling of one scattering volume. 

[0045] Figure 10a shows an example of an implementation in which multiple sampling volumes are defined in the 
25 lens LE of an eye. The front of the lens LE is positioned in upper part of Figure 10a, and the scattering apparatus is 
positioned in front of the lens, so as to measure scattering particle size at different parts of the lens in situ. A first 
approximately parallel beam B1 and a second approximately parallel beam B2 generated by first and second trans- 
mitters 10a, 10b (not shown in this Figure) are incident along the lens axis and approximately parallel to the lens surface 
respectively. Each beam may contain one or more discrete wavelengths generated by different lasers. For each beam 
30 B1 , B2 there are five receivers 12a to 1 2e which receive scattered light from respective reception volumes RV11-RV15, 
RV21-RV25 which intersect the beams B1 , B2 to define scattering volumes SV11-SV15, SV21-SV25 respectively for 
measuring particle size at different parts of the lens LE. A similar arrangement of intersecting beams B and detected 
volumes V can be used to measure particle size across a micro-vial plate. 

35 Depolarised Cross-correlation DLS 

[0046] While the translational diffusion velocity is proportional to the particle radius r, the rotational diffusion velocity 
is proportional to i 3 and is therefore more sensitive to small changes in size. Moreover, by comparison of rotational 
and translational diffusion velocities, information on particle shape can be obtained. Rotational velocity can be derived 

to from the depolarised component of scattered light. 

[0047] Figure 10b shows apparatus for depolarised cross-correlation DLS measurement. A laser 11 generates a 
beam having vertical polarisation which illuminates a sampling volume SV and the scattered light from the sampling 
volume SV is received by a receiver 12 having at least two mono-mode fibres 2 through respective parallel horizontal 
polarizing filters Fa, Fb. As in the arrangement III described with reference to Figures 1a and 1b, each fibre 2 is coupled 

45 to a detector 14, discriminator 16 and pulse former 18, and the different outputs are correlated by a cross-correlator 
20, as described with reference to Figure 10. The complementary adjustment of the polarization direction of the laser 
11 compared to the polarizing filters F removes system noise by cross-correlating the two signals. Thus, very subtle 
changes in signal strength and in the count rate of the depolarized component can be resolved. The depolarized 
component is used to calculate the time constant of rotational diffusion of the particles, which is proportional to r 3 

50 whereas the transitional diffusion time constant is proportional to r 2 , as shown by equations (3) and (5). 

Fluorescence Measurement 

[0048] In this embodiment, the apparatus as shown in Figure 1 0 is combined with the design of Figure 1 0b used with 
55 a fluorescent sample. The maximal contrast of the intensity fluctuation is observed, if the scattering volume is very 
small (e.g. 4x4x4 urn 3 ), although the fluorescent light is inelastic and incoherent. The movement of the particle 
through such a limited volume produces a quasi-stochastic signal, which can be correlated and eventually yields a 
time constant. If the particle is linked to a marker, this acquired time constant relates to kinetics or diffusion conditions 
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prevailing in the sample for the particles linked to the marker. Alternatively, if the sample is illuminated with a plurality 
of lasers of different wavelengths, autofluorescence may occur and no marker is therefore needed. 

Active Remote Sensing 

[0049] Figures 12a and 12b show two alternative embodiments for polychromatic active remote sensing of the at- 
mosphere. These embodiments represent opposite extremes of design, in which the distance between transmitter and 
receiver is respectively either very large and the scattering volume very short as in Figure 12a, or the distance is not 
very large and the scattering volume very long as in Figure 1 2b. For both embodiments, the signal processing apparatus 
is as shown in Figures 1 3a to f. 

[0050] In both embodiments, pulsed lasers 11 of respective different wavelengths are used to generate the illumi- 
nating beam. The lasers 11 are either triggered by a pulse generator 22, or is synchronised with the sampling timing 
of the correlator 20. In the embodiment shown in Figure 12b, a time-of-flight register 24 measures the time between 
the pulsing of the laser 11 and the received scattered pulse and thereby calculates the time of flight. The received 
pulse intensity from each of the detectors is input to a multi-channel integrator 27 which combines the intensity as 
received from each detector 14 and passes this to a store 28 which stores the intensity as a function of time. This 
function is passed to the processor 30 which performs evaluation of the lidar data in the conventional manner, to 
calculate scattering intensity as a function of height. 

[0051] Both the embodiments have a receiver comprising a plurality of mono-mode fibres 4, as in the arrangement 
III shown in Figures 1a and 1b, and detectors 14. Figure 13a shows the waveform of one of the pulsed lasers 11, while 
Figure 1 3b shows the received intensity signal as measured by one of the detectors 1 4. The received signal is processed 
by removing the intervals between the received signal pulses to give a concatenated signal as shown in Figure 13c. 
The concatenated signal is then digitally high-pass filtered by filter 26 to give the waveform shown in Figure 13d. The 
filtering removes the effect of fluctuation of the number of scattering particles in the scattering volume. The resultant 
signal is then cross-correlated to give the CCF shown in Figure 13e, from which the radius distribution as shown in 
Figure 13f is derived by the method described above. In the embodiment shown in Figure 12b, the autocorrelation 
function is calculated as a function of atmospheric height. 

[0052] In both Figures 12a and 12b, the optical equipment OE comprises a telescope defining the transmitted beam 
B and a separate telescope defining the receiving volume RV. 

[0053] In one example, one of the lasers is a continuous wave 1 W NdYAG laser emitting light at 532 nm is used to 
illuminate a scattering volume as given in Figure 8 at a height of 1 0000 m and a transmitter-receiver distance of 1 000 m. 

Applications 

[0054] The above embodiments, particularly the configurations of the electronics and the receiver fibres, can be used 
for many different applications, as described below. 

[0055] The apparatus shown in Figure 8 may be used both for laboratory-based studies, such as aerosol or com- 
bustion study, and for reduced gravity environment studies, such as may be carried out on spacecraft. As the apparatus 
can be miniaturized, it is particularly suitable for a spacecraft payload. Possible microgravity applications include spec- 
troscopy on complex chemical systems, fibre protein aggregations, kinetic studies and nucleation onset together with 
the study of early crystal growth. For a given payload size, the apparatus allows a higher density of devices; this is 
particularly important for a Protein Crystallization Diagnostics Facility. 

[0056] The apparatus could also be applied to most types of ground-based lidar for the determination of multiple 
scattering or sizing of particles, or performance of PCS or FCS on droplets or scattering centres, or determination of 
visco-elastic properties. 

[0057] The apparatus is particularly useful for applications in ophthalmology, where it is difficult to obtain angle- 
dependent measurements under in vivo conditions. For example, the mechanisms leading to opacification of the eye 
lens and the changes in the vitreous body as a function of age may be measured. 

[0058] Fluorescence measurements could be used to measure subtle chemical gradients in the micro-environment 
of biological systems. The apparatus for wavelength dependent particle sizing may be used without modification for 
fluorescent samples. The separation of detected fluorescent light may be accomplished using known specific aperture 
filters combined with the new arrangements shown in the third arrangement of Figures 1a and 1b. This may be used 
to achieve even lower detection thresholds, as described for example in 'Fluorescence Correlation, Single Molecule 
Detection and Large Number Screening Application in Biotechnology', Rudolf Rigler, Journal of Biotechnology 41 , pp. 
177-186 (1995)., without any sample manipulation or added markers. 

[0059] The lidars shown in Figures 1 2a and 1 2b could be installed on two satellites, with the transmitter on one, and 

the receiver on the other or the transmitter could be located on a satellite with the receiver on the ground. 

[0060] In the atmospheric sensing application, the effect of multiple scattering, the limitation on diffusion by particle 
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collision and the velocity components of particles in clouds can be differentiated. This could be used for example to 
determine the horizontal displacement of a cloud of pollutants in atmospheric boundary layers, by determining the 
vertical profile of velocity vectors. The embodiment can be used to measure properties of aerosols in the atmosphere, 
such as particle or droplet size, and tracking of pollution especially in regions of high aerosol loading. The determination 
of particle sizes in the troposphere is particularly important, as ozone depletion is thought to be linked to the nature of 
particles in tropospheric clouds ('Unexpected Low Ozone Concentration in Midlatitude Tropospheric Ice Clouds: A 
Case Study', J. Reichardt et. al., Geophysical Research Letters, Vol. 23, No. 15, pp. 1929-1932, 15 July 1996). 
[0061 ] The active remote sensing embodiment of the present invention may be used to determine in situXhe pressure 
and temperature dependency of condensation of atmospheric water vapour. The measurement of smallest particle 
sizes as a function of atmospheric height may be used to estimate water vapour cross-sections, for which the published 
data are in considerable disagreement. The application to atmospheric measurements allows the measurement of the 
wavelength dependency of scattering by ice-crystals or droplets in optically thick layers. Dutch patent no. 1001369 
describes how DLS or PCS test results may be applied to data derived from low concentrations of particles suspended 
in air as aerosols. 

[0062] Embodiments of the present invention allow a correction factor to be determined between multiple scattering, 
the scattering coefficient and the diameter adjustment of the field of view. The correction factor is determined by cal- 
culating the separation between the theoretical curve and the measured curve for the sample under observation, as 
shown for example in Figure 7. 

[0063] For example, if multiple scattering occurs, the derived ACF will not be plotted as a function of intercept and 
sampling time structure of the correlator, but as a function of the square root of time. This is called diffusive wave 
spectroscopy, as described in Chapter 16, 'Diffusing Wave Spectroscopy' by David A. Weitz and David J. Pine of the 
book 'Dynamic Light Scattering - The Method and Some Applications' edited by Wyn Brown, Clarendon Press 1993, 
and could be applied to the above embodiments in real time. 



Claims 

1. A dynamic light scattering method, comprising: 

illuminating a fluid scattering volume (SV) with temporally and spatially coherent radiation (B), and 
detecting the fluctuation of radiation scattered by said fluid volume (SV), 

characterised in that the fluid scattering volume is illuminated with said radiation at a plurality of discrete 
wavelengths, and 

in that the detecting step comprises detecting the intensity fluctuation of each of said wavelengths of said 
scattered radiation. 

2. A method as claimed in claim 1 , further comprising correlating the detected intensity at each of said wavelengths 
as a function of time so as to derive a property of scattering centres in said fluid scattering volume (SV). 

3. A method as claimed in claim 1 or claim 2, wherein the fluid scattering volume (SV) is simultaneously illuminated 
with each of said plurality of discrete wavelengths. 

4. A method as claimed in any preceding claim, wherein each of said wavelengths is detected at approximately the 
same scattering angle. 

5. A method as claimed in any preceding claim, wherein the scattered radiation is received by an optical arrangement 
including a mode-selective transmission medium (2) prior to detection. 

6. A method as claimed in claim 5, wherein the scattered radiation is received by a plurality of mono-mode optical 
fibres (2) each coupled to a respective receiver (14) for detecting the radiation received by that fibre (2). 

7. A method as claimed in claim 6, wherein the intensity fluctuation of the scattered radiation of one or more of said 
wavelengths is received by at least one pair of said receivers (14) corresponding to each of said one or more 
wavelengths and the outputs of each said pair are cross-correlated. 

8. A method as claimed in any preceding claim, wherein at least one further fluid scattering volume (SV) is illuminated 
with said radiation (B), and the intensity fluctuation of said radiation scattered from said at least one further fluid 
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scattering volume (SV) is detected. 

9. Dynamic light scattering apparatus, comprising: 

a radiation source (11) arranged to illuminate a fluid scattering volume (SV) with temporally and spatially co- 
herent radiation (B), and 

a detector (14) arranged to detect the fluctuation of said radiation scattered by said fluid volume (SV); 

characterised in that said the radiation source (11) is arranged to generate said radiation (B) at a plurality of 
discrete wavelengths, and 

in that the detector (14) is arranged to detect the intensity fluctuation of each of said wavelengths of said 
scattered radiation. 

10. Dynamic light scattering apparatus, comprising: 

a radiation source (11) arranged to illuminate a fluid scattering volume (SV) with temporally and spatially co- 
herent radiation (B), and 

a detector (14) arranged to detect the fluctuation of said radiation scattered by said fluid volume (SV); char- 
acterised in that 

the detector (14) comprises a plurality of mono-mode fibres (2) each arranged to detect said radiation at 
approximately the same scattering angle and each coupled to a respective radiation detecting element. 

11. Apparatus as claimed in claim 10, wherein the radiation source includes one or more lasers (11) each coupled to 
a respective one or more mono-mode fibre (4). 

12. Apparatus as claimed in claim 10 or claim 11, wherein the radiation source (10) transmits said radiation through, 
and the detector (14) receives said scattered radiation through, respective different optical elements (OE). 

13. Apparatus as claimed in any one of claims 10 to 12, including a cross-correlator (20) for cross-correlating as a 
function of time the outputs of at least two of said radiation detecting elements (14) receiving said scattered radiation 
at the same wavelength. 

14. Apparatus as claimed in any one of claims 10 to 13, wherein said illuminating radiation (B) is polarized, and the 
detector (14) includes polarizing means (Fa, Fb) so as to detect a depolarized component of said scattered radi- 
ation. 



15. A method of detecting the velocity of scattering centres in a fluid scattering volume (SV) by means of dynamic light 
scattering, comprising: 

illuminating the fluid scattering volume (SV) with temporally coherent radiation (B) of a discrete wavelength, 
detecting the fluctuation of radiation scattered by said fluid volume (SV), 

correlating the detected intensity as a function of time to generate a correlation function (ACF-33); 
determining one or more oscillation frequencies of components of said correlation function (ACF-33); and 
deriving velocity components of scattering particles in said fluid volume (SV) from said one or more oscillation 
frequencies. 

16. A method as claimed in claim 15, wherein the velocity components are derived by means of the relation 



where v e is the velocity component is the direction of detection, Af is the oscillation frequency, f is the fre- 
quency of the radiation, n is the refractive index of the fluid volume, and c is the speed of light. 

17. A method of measurement of properties of a scattering volume (SV) of the earth's atmosphere by dynamic light 
scattering, comprising: 

illuminating said scattering volume (SV) of the atmosphere with pulses of temporally and spatially coherent 
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radiation , 

detecting the fluctuation of said radiation pulses scattered by said fluid volume (SV), and 

correlating said detected fluctuation as a function of time so as to derive said properties of the scattering 

volume (SV). 

18. A method as claimed in claim 17, further comprising: 

concatenating the detected pulses to generate a fluctuation function; wherein said correlating step comprises 
correlating said fluctuation function. 

1 9. A method as claimed in claim 1 8, further comprising high-pass filtering said fluctuation function prior to correlation. 

20. Apparatus for measurement of properties of a scattering volume (SV) of the earth's atmosphere by dynamic light 
scattering, comprising: 

a radiation source (11) for illuminating said scattering volume (SV) of the atmosphere with pulses of temporally 
and spatially coherent radiation, 

a detector (14) for detecting the fluctuation of said radiation pulses scattered by said fluid volume (SV), and 
a correlator (20) for correlating said detected fluctuation as a function of time so as to derive said properties 
of the scattering volume (SV). 

21. A fluorescence correlation spectroscopy method, comprising: 

illuminating a fluid volume including fluorescent particles with radiation, and 
detecting the intensity fluctuation of fluorescence emitted by said fluid volume, 

characterised in that the fluid volume is illuminated with said radiation at a plurality of discrete wavelengths, 

and 

in that the detecting step comprises detecting the intensity fluctuation at each of said wavelengths of said 
fluorescence. 

22. A method as claimed in claim 21 , further comprising correlating the detected intensity at each of said wavelengths 
as a function of time so as to derive a property of fluorescent centres in said fluid volume. 

23. A method as claimed in claim 21 or claim 22, wherein said fluorescent particles are substantially free of fluorescent 
markers. 
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